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ABSTRACT 

We present a deep near-infrared wide-angle photometric analysis of the structure of 
the inner Galactic bar and central disk. The presence of a triaxial structure at the 
centre of the Galaxy is confirmed, consistent with a bar inclined at 22±5.5° from the 
Sun-Galactic centre line, extending to about 2.5 kpc from the Galactic centre and 
with a rather small axis ratio. A feature at £=—9.8° not aligned with this triaxiality 
suggests the existence of a second structure in the inner Galaxy, a double triaxiality 
or an inner ring. We argue that this is likely to be the signature of the end of the 
Galactic bar, at about 2.5-3 kpc, which is circumscribed by an inner pseudo-ring. No 
thick dust lane preceding the bar is detected, and a hole in the disc's dust distribution 
inside the bar radius is inferred. 

Key words: Galaxy: structure - Galaxy: stellar content - Galaxy: bulge - galaxy: 
bar - infrared: stars - dust, extinction. 



1 INTRODUCTION 

The central region of the Galaxy is our closest opportu- 
nity to quantify the structure and formation of a dense and 
massive complex stellar-dynamical system. A detailed star- 
\ by-star analysis is possible in the Milky Way bulge, in con- 
trast to distant galaxies where only the integrated light is 
observed. But our point of view from within the Galactic 
disc has its drawbacks. The high interstellar extinction, the 
crowding and the confusion between foreground disc stars 
and bulge sources, make studies of the inner Galactic re- 
gions difficult. 

There is substantial evidence for the presence of a tri- 
axial structure in the inner Galaxy. The presence of a bar 
was first suggested from the obs ervations of non-circula r mo- 
tions in the gas kinematics bv lde Vaucouleursl l|l964j) . The 
first direct evidence came from detection of an asymme- 
try in the infrared luminosity distributi on jBlitz fc Spere ell 
Il99llh confirmed by the COBE map JPwek et al.lll995l) . 
Star counts confirmed the presence of an asymme- 



iNakada et alJll99ll.lWeinbere 


f992), near-infrared star sur- 


vevs JHammerslev et al.ll994l. 


LoDez-Corredoira et alJl997l 



of red clump stars dStanek et alJ Tf994t) . The new near- 
infrared surveys DENIS and 2MASS, with arcsecond spa- 
tial resolution and sufficient sensitivity to detect typi- 
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cal red giant bulge sources, complemented by the mid- 
infrared IS O survey ISOGAL, are beginning to provide new 
const raints llLopez-Corredoira et al.ll200ll . ICole fc Weinberg 



l200l Ivan Loon et al.ll200St . The high optical depth mea- 
sured by the microlensing surveys also points towards a non- 
axisymmetric model, but the exac t constraints remain con- 
troversial (e.g. iBinnev et al. 2000). 

Only a few dynamical constraints are available so 
far, essentially from stellar kinemat ics from SiO m asers 
JPeguchi et alJ l200Ct> . OH/IR stars (ISevensterlll999lh and 
from observations through a few low-extinction windows 
(Hafn er et alJl200Cf ) . Early dynamical m odels of gas HI, CO 
and CS emissions feinnev et alJll99ll) and infrared lumi- 
nosity distribution models jDweket^l]|l995l iFreudenreicbl 
1998) were derived independently, while more recently new 
models are being developed to incorporate both gas and lu- 
minosity c onstraints, using in particular N-body simulations 
jFuxfll999lh 

All those studies agree that the inner Galaxy at low lat- 
itudes is asymmetric, with the bulge/bar major axis at pos- 
itive galactic longitudes. The exact orientation, length and 
axis ratio of the inner asym metric structure are still poorly 
constrained <Gerhardl2 001h The relative importances of the 
bulge and the inner disk in the central regions, the stel- 
lar population content, a nd the formation and evolutionary 
history remain unknown JWvse et alJll997lh Discussion con- 
tinues on the need f or a distinctio n between a central bulge 
and a long thin bar iKuiikenll996t) . Indeed, along the line of 
sight towards the central regions of the Galaxy one crosses 
the galactic disc, spiral arms, the molecular ring, the bar 
and/or bulge, and maybe other complex structures such as 
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a hole in the disc distribution, or a possible stellar ring or 
double bar, as observed in external galaxies. To improve our 
understanding of this crucial region new high sensitivity and 
high spatial resolution multi-wavelength studies are needed. 

In this paper we present new deep near-infrared pho- 
tometry of the inner Galaxy close to the Galactic plane, 
taken with the Cambridge InfraRed Survey Instrument 
(CIRSI). Those observations provide large samples of 
bulge/bar red clump stars, allowing extinction mapping and 
reliable 3-D spatial analyses in fields at both positive and 
negative galactic longitudes in the Galactic plane. 



2 CIRSI PHOTOMETRY OF THE INNER 
GALAXY 

The requirements for an observational study of inner Galac- 
tic stru cture are well-established , following earlier stud- 
ies re.g.lUnavane fc OilmorellT^ Hhata, fc Gilmorelll 995al. 
llbata fc Gilmorelll995bri . Near infrared observations are op- 
timal, to identify the dominant intrinsically-bright red giant 
stellar population, while minimising the effects of interstel- 
lar extinction (longer wavelength mid-infrared studies, such 
as those for ISOGAL, are optimised for study of the latest M 
giants). Observations of the inner Galaxy preferably in the 
Galactic plane, and certainly at low Galactic latitudes, are 
desirable to minimise model-dependent corrections for the 
as-yet poorly known scale heights of the inner-galactic stel- 
lar populations. To study the expected asymmetry caused by 
the Galactic bar, observations at both positive and negative 
longitudes are essential, with, in so far as is possible given 
the ext inction distribution, symm etric longitude fields pre- 
ferred. lUnavane fc Gilmord Jl99S) have shown that struc- 
tural studies at absolute galactic longitudes greater than 4 
degrees are needed to constrain bar models to a useful ex- 
tent. 

Given these requirements, we selected for further de- 
tailed study fields observed by DENIS and ISOGAL at 
galactic latitudes —0.3° < b < 0.3° and galactic longitudes 
I = ±5 — 6°, where the effect of any bar-like feature should 
be strong, and at £=±9 — 10° where the bar may or may 
not be present. The field C32 (£ = 0° ,b = 1°), which was 
used as an ISOGAL calibration field be cause of its low and 
uniform extinction llOmont et al]ll999ft . was also selected 
for calibration purposes here, and in addition to provide the 
minor-axis reference to set the zero-point of the expected 
asymmetry. 

Near-infrared observations in the J(1.25/im), 
H(1.65/im) and K s (2.15/im) bands were obtained with 
the Cambridge Infrared S urvey Instrument jBeckett et alJ 
Il997l. iMackav etai]|200rJl on the du Pont 2.5m telescope 
at Las Campanas Observatory. CIRSI is a mosaic imager 
consisting of four Rockwell IK x IK detectors. The pixel 
scale is 0.2 arcsec pixel -1 . The gaps between the detectors 
being comparable to the detector size, four dither sets are 
used to create a fully-observed mosaic image, leading to a 
basic field of view of about 13 x 13 arcmin 2 . 

The observations carried out for this survey are sum- 
marised in tabled 



2.1 Data reduction 

The data reduction was carried out using an exten- 
sively updated version of the InfraRed Da ta Reduction 
(IRDR ) software package, first developed bv ISabbev et alJ 
(2001). A summary of the full process is given here. 
The updated version of IRDR with its documentation 
was developed by one of us (CB) and is available at 
|http:/ /www.ast. cam.ac.uk /~optics /cirsi/software 

First each image is corrected for non-linearity, as the 
array detectors used are non-linear in their response to flux. 
To calibrate the non-linearity, domeflats were taken with 
different exposure times, with a short exposure observation 
between each as reference. The expected flux for an exposure 
time is computed from the reference exposure flux level. The 
ratio of this flux to the measured flux is observed to increase 
quadratically with the flux level. A linear quadratic regres- 
sion is computed up to the saturation level (about 40,000 
ADU), setting the first coefficient to unity to constrain no 
correction for a zero signal level. This relation is then used 
to correct the non-linearity on all images. 

The detectors generate an internal (thermal noise) sig- 
nal even when there is no external signal. This 'dark' signal 
must then be subtracted 1 . This is a straightforward process, 
using exposures of the same length as the science exposures, 
but with the shutter closed. 

The data are then flatfield corrected. This process uses 
the difference image derived by subtracting domeflats ob- 
tained with the illumination lamp turned on from sub- 
sequent domeflats obtained with the illumination lamps 
turned off. This difference image is normalized to the sensi- 
tivity of the first detector. These flatfields are also used to 
detect bad pixels and to create weight maps, which are used 
to create a signal-to-noise value for each pixel observation 
for use during the coaddition of all the individual 'dither' 
exposures into a final single combined image. 

The sky is subtracted in two passes. The basic obser- 
vational 'unit' is a set of typically five 20-second repeated 
exposures ('loops') on a single field. Nine of these sets of ex- 
posures are obtained, with small telescope pointing offsets 
('dither') between each set. The five loops are combined, us- 
ing sigma clipping, into a single frame called 'dither frame'. 
A first-pass sky image is derived by median-combining the 
nearest dither frames. This process of course leaves the real 
sources in the frame. The nine dither frames are then off- 
set and a preliminary combination is made. After this first 
dither frame coaddition, obje ct masks are produced using 
SExtractor source extraction jBertin fc Arnout j[l993) and 
used to calculate a mask which can exclude all detected 
sources from the raw data frames. Object-masked frames 
are used to make a second pass sky subtraction on each ba- 
sic loop image. 

The spatial offsets generated by the dithering between 
the dither frames are computed by cross-correlating object 
pixels detected by SExtractor. The nine individual dither 



1 During the observations made at the end of September and 
during October 2001, a nitrogen leak from the dewar, with conse- 
quent temperature drifts, affected the stability of the dark count 
and the linearity of the chip 3 detector so badly that those data 
had to be discarded. 
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Table 1. Log of the CIRSI observations. Positions of the centre of the mosaic are given both in galactic and equatorial J2000 coordinates. 
The exposure mode is given by the number of dither frames x the number of sub-integrations(loops) with their individual exposure time. 
( 1 ) after the exposure time indicates that chip 3 has been discarded, see footnote 1 . The seeing conditions are indicated by the PSF 
FWHM with its variation range during the mosaic observation. The last column indicates the effective useful magnitude limit. 



Field (i,b) 


RA (J2000) 


Filter 


Date 


Exposure 


seeing 


mag limit 












\ ) 


(mag) 


5NN (-5.73,-0.22) 


17:32:08.0 


J 


2000-08-17 


9 x 5*20s 


0.6-0.64 


21.0 




-33:54:00.0 


H 


2001-09-06 


9 x 5*20s (!) 


0.72-0.84 


18.3 






K 5 


2001-04-10 


9 x 3*20s 


0.68-0.7 


18.3 


5NP (-5.74,+0.19) 


17:30:30.0 


J 


2000-08-17 


9 x 5*20s 


0.66-0.72 


20.9 




-33:41:00.0 


11 


2000-08-16 


9 x 5*20s 


0.66-0.82 


19.5 






-tVs 




y X Z OUa 


nan ao. 


la 


5PN (+5.67,-0.28) 


17:59:35.0 


J 


2000-08-17 


9 x 5*20s 


0.8-0.96 


20.2 




-24:12:00.0 


11 


2000-08-16 


9 x 5*20s 


0.6-0.76 


20.2 






|<- 

*Vs 


onm rizi in 

aUU1-U%- jlu 


y x o zus 


n a_a n q 

U.DD-U.y 


lo.z 


5PP (+5. 76, +0.23) 


17:57:53.0 


J 


2001-09-04 


9 x 5*20s 


1.04-1.24 


19.6 




-23:52:00.0 


11 


2000-08-16 


9 x 5*20s 


0.76-0.8 


19.3 






K s 


2001-04-15 


9 x 2*30s 


0.68-0.84 


18.4 


9N (-9.80.+0.05) 


17:19:51.6 


J 


2001-09-07 


9 x 5*20s ( x ) 


1.04-1.2 


19.7 




-37:07:21.9 


H 


2001-09-08 


9 x 5*20s (!) 


0.8-0.84 


18.9 






K s 


2001-09-03 


9 x 5*20s 


0.8-0.86 


18.6 


9P (+9.55,-0.09) 


18:07:09.6 


J 


2001-09-04 


9 x 5*20s 


1.4-1.7 


19.4 




-20:43:28.5 


H 


2001-09-06 


9 x 5*20s (!) 


0.86-0.96 


18 






K s 


2001-09-03 


9 x 5*20s 


1.06-1.4 


17.9 


C32 (+0.00,+1.00) 


17:41:45.0 


J 


2001-09-07 


9 x 5*20s (!) 


0.94-1.1 


18.1 




-28:25:00.0 


H 


2001-09-08 (3/4) 


9 x 5*20s I 1 ) 


1.06-1.3 


17.7 






H 


2001-10-02 (1/4) 


9 x 5*20s I 1 ) 


0.86 








K s 


2001-09-05 


9 x 5*20s I 1 ) 


0.9-1.1 


16.6 



frames are then coadded using a weighted bi-linear interpo- 
lation, excluding bad pixels. 

Finally, the astrometry is calibrated by correlating the 
SExtractor's object catalogue with the 2MASS catalogue. 

In practice, the image PSF is seen to vary significantly 
during the typically 15 minutes which are required to create 
a full set of dithered observations. Consequently, to max- 
imise the photometric quality of the data, no attempt was 
made to create single full mosaic images. 

2.2 Point source extraction and photometry 

The borders of the images, which do not have all the dither 
observations and so have highly variable signal-to-noise ra- 
tios, are removed before beginning the source extraction pro- 
cedure. 

PSF-fitting photometry was carried out with the IRAF 
DAOPHOT package. A first list of detections was created 
by the DAOFIND procedure with a 5— a threshold. An ini- 
tial estimate of the photometry is then provided by the 
PHOT procedure, adopting as aperture the Full Width at 
Half Maximum (FWHM) value calculated from the image 
Point Spread Function (PSF) (see Table QJ. A list of bright 
and 'almost' isolated stars, preselected by the PSTSELECT 
routine, was interactively confirmed. Those stars are used 
by the PSF task to compute a PSF model, using a Gaussian 
profile and a lookup table quadratically varying with the po- 
sition in the image. The PSF model is then refined by itera- 



tively subtracting the neighbours of the stars used to define 
the PSF, before computing the final iterated PSF model. A 
single PSF FWHM value is computed for each 'dither set', 
composed of the 4 detectors images, but a separate spatially- 
varying PSF model is calculated for each detector image. 
Finally the ALLSTAR procedure fits the model to all the 
stars detected. A second-pass detection is performed on the 
residual image created by subtracting all sources detected 
in the first-pass. This second-pass uses a higher threshold of 
8— cr, to avoid detecting subtraction residuals. 

DAOPHOT provides two measures of the reliability of 
a detection and of its photometric accuracy: the goodness of 
fit, x 2 , an d a measure of the image sharpness, sharp, indicat- 
ing image blemishes and resolved objects. Detections with 
X 2 > 3 and | sharp] > 2 have been eliminated. A more severe 
selection of stars with reliable photometry is described be- 
low. Detections around highly saturated stars are manually 
deleted. 

Standard stars from lPersson et ail dl998t) were observed 
to derive the magnitude zero-point for each night. They were 
reduced as above, and analysed to derive zero-point photom- 
etry using an aperture photometry radius of 20 pixels, equiv- 
alent to a diameter of 8 arcsec. Those observations show that 
a large number of the fields of this survey where observed 
during non-photometric nights. We thus require an exter- 
nal calibration in these fields, which we derive from 2MASS. 
The H and K s filters of CIRSI and 2MASS are identical, 
ensuring a straightforward relative calibration. However the 
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Figure 1. Estimation of the photometric errors in the K s band 
from repeat observations with a different central pointing in field 
5PN. 
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Figure 2. Position of all stars detected in field 5NP in J H and 
K s with x 2 < 1-5 and \sharp\ < 0.5. 



J band filter of 2MASS is more extended into the atmo- 
sphcri c water ab sorption features at around 1.1 and 1.4 /im 
JCarpenterll200lD . leading to a photometric zero-point shift 
of typically J^cirsi) - J(2MASS) ~ -0.04 mag. To ensure 
that all our photometry is on a homogeneous photometric 
system, we then calibrated all our fields to the 2MASS cat- 
alogue zeropoints. The calibration was done on each dither 
set, the atmospheric conditions changing on this time scale 
during non-photometric nights. By comparing this 2MASS 
calibration to the standard star calibrations for good nights, 
the zero-point calibration accuracy is better than 5%. 

To further test the accuracy of our data reduction and 
source extraction, we obtained a fifth full data set for field 
5PN in the K 3 band. This field was pointed at the centre 
of the standard data sets for this field, and so provides an 
independent data set which overlaps all the other dither sets. 
This independent set of photometry was then analysed to 
determine the difference in the photometry for this field as 
a function of magnitude. The resulting dispersion, presented 
in figure should be conservative as the overlapping areas 
are at the corner of the images where the PSF fitting is the 
least robust. On the other hand those dispersions could be 
larger during non-photometric nights. 

Combining the imaging data for the different filter pass- 
bands is straightforward. Each single-colour image is divided 
into several small areas. In each area for each colour, stars 
with very reliable photometry and astrometry {\ 2 < 1.5 and 
\sharp\ < 0.5) are cross-correlated between different pass- 
band images, and matched with an allowed matching radius 
of 1.3 arcseconds. The mean shifts between the astrometric 
systems of the two passbands are then computed and ap- 
plied to all stars. A final cross-identification is then defined, 
using these adjusted astrometric solutions, with a smaller 
matching radius of 0.5 arcseconds. The JHK S catalogue is 
then derived from the JH and JK S catalogues, using the HK S 
cross-match to reject false and multiple cross-identifications. 

Given the complexity of the source distributions in our 
fields, a completeness estimation for number counts analyses 
would be unreliable. Indeed figure [5] illustrates the various 
instrumental, atmospheric and astrophysical effects on the 
number of stars detected. The gaps between the detectors 
are apparent due to our removal of the dithered image bor- 



ders. The photometric conditions influence the overall com- 
pleteness. Variations of the number counts between dither 
sets are due to changes in the atmospheric conditions and in 
particular in the seeing. Large holes are due to the presence 
of saturated stars. Variations of the number counts within 
chip images are also visible due to variations of the extinc- 
tion. For example a dark cloud is clearly present in the top- 
left corner of figure [5] 



3 RED CLUMP GIANTS AS DISTANCE 

INDICATORS: DERIVING THE DISTANCES 

The main stellar populations observed in this survey are il- 
lustrated in figure|3] which is the colour-magnitude diagram 
(CMD) of the minor axis low-extinction standard field C32. 
In the blue part of the diagram lie foreground main-sequence 
disc stars. In field C32, since the extinction is low and be- 
cause the central inner Galactic bulge clearly is spatially 
highly-concentrated with a small scale length, the bulge red 
giant branch (RGB) stars form a well-defined, almost single- 
distance RGB feature with a distinct red clump. The fore- 
ground disc red clump giants form a third sequence on the 
CMDs, roughly parallel to the disk main sequence stars, be- 
coming fainter and redder as their distance and total extinc- 
tion increase. The narrow intrinsic luminosity distribution of 
the red clump stars, apparent in this figure, illustrates their 
utility as good distance indicators. 

Given the high and variable extinction apparent in our 
fields, reddening independent magnitudes are a critical re- 
quirement. Such parameters can be derived using any two 
colours, given an extinction law, for example: 

K SJ _ Ks =K s - Ak ' (J-K s ) (1) 
Aj — Ak s 

In the near-infrared, the extincti on curve is well de- 
fined by a single power law fu nction (e.g. lCardelli et alll98St 
lMathill990l . lHe et all 199^ : A x /Aj = ( A/A j) "". As a con- 
sequen ce, and unlike optical studies such as IStanek et alJ 
(1994), where extinction law variations remain a source of 
possible systematic uncertainty, our reddening independent 
magnitudes should not be significantly affected by variations 
of ISM properties along the different lines of sight, although 
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Figure 3. Colour-magnitude diagram of the bulge field C32 (£ = 
0°, b = 1°), using all stars detected in J,H and K s with \ 2 < 1-5 
and \sharp\ < 0.5. The dotted arrow represents the movement of 
a star on this diagram due to an absorption of Av=5 mag. 



estimates of Ay w ould. In the foll owing we adopt the coef- 
ficient a = 1.73 of iHe et alJ (119951) . leading to Ak, IEj-k, 
= 1X64, which agrees with the values reviewed by [Mathis 
l|l99(t, . 

In a colour magnitude diagram made with a reddening- 
independent magnitude, the effect of distance and extinction 
are disentangled: a star of a given spectral type moves ver- 
tically with distance, horizontally with extinction. Figure 0] 
illustrates, for each of our seven observed fields, how the lo- 
cation of a red clump star in the (J — K s , K s j-k s ) diagram 
can be translated into an (Ay, distance) estimate. These dia- 
grams therefore allow us to quantify the distribution of both 
stellar density and of the extinction along each line of sight. 
For example, in field 9P a clear flattening of the disc red 
clump distribution at K sJ _ Ks « 10.5 indicates the presence 
of a specific source of extinction located at about 3.5 kpc 
from the Sun, which increases Ay by about 7.5 mag in only 
1.5 kpc. This feature is consistent with the steep increase 
in abs orption at 3.4-4 kpc seen at £=12.9° bv lBertelli et alJ 
(1995), that they associated with the molecular ring. In all 
our galactic plane fields, the bulk of bulge red clump stars 
become visible at Av ~ 10 mag, and apparently suffer a 
range of internal extinction of more than AAv ~ 5 mag. 

A clear difference between the red clump distance dis- 
tributions is present between the different lines of sight, in- 
dicating that we have resolved the spatial structure of the 
inner Galaxy. 

To quantify this, we determine the reddening- 
independent magnitude of the red clump, mnc, by a 
non-linear least-squares fit of functi on |21 to the his- 
togram of red giant st ars (e.g. IStanek fc Garnavichl Il99ct 
ISaUria fc GimrTflEool): 

ati ^ i i i 2 . Nrc r {iriRc - m) 2 

N (m) = a + bin + cm H -= exp[— ^ — J (2) 

The Gaussian term represents a fit to the bulge/bar red 
clump. The first three terms describe a fit to the background 
distribution of non-bulge-clump red giant stars, which in our 
case contains not only the other bulge giants, but also the 
disc red clump stars that fall into the bulge red clump. The 



foreground dwarf star sequence is eliminated on the CMDs 
by the criterion K s > 4(J — K s ) + b, with b determined for 
each CMD to account for the different extinction. As dis- 
tance and extinction grows, the mixing between faint dwarfs 
and giants increases. Only stars with K SJ _ Ka < 14 mag will 
then be used, which corresponds to a limit in red clump star 
distance of about 15 kpc. These selection criteria ensure that 
all red clump stars at a distance of 9 kpc with photometric 
error smaller than 3 a and extinction smaller than ^4^ = 15 
are included in the analysis. 

For given values of the red clump absolute magnitude 
and colour, a reddening independent magnitude provides a 
direct estimate of the distance modulus /i of red clump stars 
(equation |3J , and so a direct measure of the line-of-sight 
variation in Galactic structure. 

To check our calibration of the red clump colours, de- 
scribed below, and any effects of photometric incomplete- 
ness, we have estimated the distance modulus of the red 
clump from all three available reddening-independent mag- 
nitudes: K 3j _ Ks , Hh-k s and Jj-h- 

mc = K sJ _ Kb + A K /Ej- K (J - K s )o - M 0Kg 
= H H -k s + A„/E„- K (H - K s ) - M 0H 
= Jj-h + Aj/Ej- H (J-H) -M j (3) 

These results are shown in figure 

We now require an upper limit to the reddening- 
independent magnitudes that can be analysed reliably in 
the fit to equation |3 For this, we derive a rough estimate 
of the relative completeness of each catalogue for each field 
(Jc ,Hc ,Kc) by determining when the number counts as a 
function of magnitude stop increasing. The different atmo- 
spheric, instrumental and astrophysical effects influencing 
this relative completeness have been described in section 2. 
The relative completeness of each reddening independent 
magnitude is then derived by combining these complete- 
ness estimates for each magnitude and field, using the mean 
colour of the giants ((J — K)m) and adding a margin of 0.5 
mag to take into account the spread in the mean colour: 

max{Kj- K ) = min(Jc — {J — K) M , Kc) 

—Ak/Ej-k * (J — K)m + margin 

The corresponding completeness-induced limits to our 
distance determinations down each line of sight are indi- 
cated for each passband and field in figure |H| by vertical 
lines. A maximum value of 16.1 is also apply due to the 
limit K sJ _ Ks < 14 mag set before for the selection of giant 
stars. 

The model fits to the spatial distribution of bulge/bar 
red clump stars down each line of sight are presented in fig- 
urc|S] In general, the fits using the three different reddening 
independent magnitudes are all consistent within about 0.1 
mag, expect for field 5NP where the dispersion is 0.2 mag. 
However no fit was possible using Jj-h for fields 9P and 
9N. 

The exception is the two fields at positive longitude 
I = 5°: no fit converged for any of the distance modulus 
estimators for fields 5PP and 5PN. The large spread in the 
red clump colour and magnitudes observed in those fields 
(figure^J cannot be explained as an artefact of completeness 
nor by large photometric errors. The interpretation of these 
results is discussed further in the next section. 
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Field 5PP Field 5NP 




Field 5PN Field 5NN 





Field C32 




Figure 4. Colour magnitude diagrams using reddening independent magnitudes of all stars observed in J and K s . In these figures, 
increasing extinction moves a star horizontally to the right, while increasing distance moves it vertically downwards. Locations of red 
clump stars on this diagram as a function of absorption and distance are given by the top and right hand axes. 
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Figure 5. The distance modulus distributions down each line of sight, computed for red clump stars using reddening-independent 
magnitudes. The vertical lines show estimates of the photometric completeness limits. Fits to equation [5] are overlaid, where available. 
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Table 2. This table presents both the mean and the line-of-sight 
dispersion in distances of bulge red clump stars for four different 
galactic longitudes. <r m is the Gaussian dispersion of the fits of 
cquation|5]to the photometric data of figure[I] which are shown in 
figure [S] <7£) is our estimate of the equivalent spread in distance, 
after deconvolution with the intrinsic red clump luminosity dis- 
persion and photometric errors. 



1 


Urc (mag) 


D (kpc) 


<trC (mag) 


a D (kpc) 


-9.8° 


14.90±0.04 


9.5±0.2 


0.25±0.03 


0.6 


-5.7° 


14.97±0.04 


9.9±0.2 


0.40±0.01 


1.5 


0.0° 


14.51±0.03 


8 


0.30±0.05 


0.8 


+9.6° 


13.63±0.04 


5.3±0.1 


0.65±0.1 


1.5 



In order to convert our reddening-independent magni- 
tudes into true distance moduli, and hence into true spatial 
density distributions, we need to calibrate the near-IR in- 
trinsic luminosity of the red clump stars. We used the true 
intrinsic colours of red clump star s predicted by the P adova 
isochrones in the 2MASS system feonatto et al-l feoCM) for a 
10 Gyr old population of solar metallicity: (J — i\" s )o=0.68 
and (H — K s )o=0.07 . We then calibrate the absolute magni- 
tude of the red clump Mok„ using our data for field C32. The 
colour-magnitude diagram for field C32, shown in figure 
is dominated by the bulge and located at 0° in longitude, its 
red clump stars should therefore be at the distance of the 
Galactic centre. If we assume the absolute r ed clump dis - 
tance of M K = -1.61 ± 0.03 mag, derived bv lAlvesl (120001) 
for the Hipparcos red clump, we obtain a distance for the 
Galactic center of Dec = 7.6 ± 0.15 kpc. Comparing the 
Padova iso chrones provi d ed in the Johnson-Cousins-Glass 
system by iGirardi et all <2002l) and the 2MASS on es, we 
find th at no correction need to be applied between the lAlve'sl 
(2000) K pho tometry and our K 3 one for red clump stars. 
According to ISalaris fc Girardil a2002r>. a small population 
correction should be applied to the lAlvea (|2000T) calibration 
for bulge stars. Using their population correction for Baade's 
window with solar metallicity we derive Mk = —1.68 mag 
and Dec = 7.8 kpc, while with enhanced a-elements it leads 
to Mk = —1.72 mag and Dqc = 8.0 kpc. However if we as- 
sume a-enhancement, as we will see later on, we also have 
to change the assumed (J — K s )o, which leads us back to 
Dgc = 7.7 ±0.15 kpc. Those values are consist ent with the 
latest Galactic centre distance estimates (e.g. iReidl Il993t 
iMcNamara et al.l l200cl. lEisenhauer et all 120031 which give 
Dgc = 8 ± 0.5 kpc. Considering the different uncertainties 
in the absolute magnitude of the red clump, we decided to 
calibrate this latter to M^f = —1.72 mag, assuming a dis- 
tance for the Galactic centre of 8 kpc. 

The results of our different fits to equation [5] have been 
averaged and are presented in table [5] The dispersion of 
the Gaussian fitted to equation |5] orc , is the convolution 
of the true line of sight dispersion in distance of the red 
clump stars (an), the intrinsic dispersion of t he red clump 
luminosity (<ro) and photometric errors (cr e ). lAlvesI <|2000T) 
estimates <to to be around 0.15-0.2 mag. The photometric 
errors are between 0.05 and 0.1 mag. An estimate of the 
dispersion due to the variation in distance modulus is then 
indicated in table [5] derived from the simple deconvolution 
on = \Jo\ c al. 



We studied the robustness of our results against differ- 
ent choices of isochrones. A different age changes the ab- 
solute magnitude of the red clump but not its colour, so 
as we fix M% c from our data the choice of age does not 
affect our results. However metallicity and enrichment in a- 
elements does affect the colours. A change in [Fe/H] of 0.4 
dex changes J—K 3 by about 0.08 mag, the higher the metal- 
licity the redder the colour. Isochrones for Q-enhanced stars 
lead to J — K s bluer by about 0.1 mag. However as all is cal- 
ibrated on field C32, the resulting distance estimates do not 
change by more than 0.3 kpc. However we note that the three 
different distance modulus estimates for field C32 (bottom 
of figure [^J agree much better with a-enhanced isochrones 
(within 0.01 mag) than with the basic set (where the disper- 
sion is 0.1 mag). On the other hand, for all the other fields 
the a-enhanced isochrones lead to a larger dispersion. If we 
assume that we are not probing one but two different stellar 
populations, which would correspond to the bulge for field 
C32 and to a distinct Galactic bar for the other fields, such a 
difference would b e expected. Chosing the ct-enhanced colors 
for th e bulge field (iMcWilliam fc Richll99llMatteucci et alJ 
1999) and the basic set for all the others would lead to an 
increase in distance modulus in all the 'bar' fields of about 
0.1 mag. I f we also consider an ag e difference, using an age 
of 6 Gyr JCole fc Weinberell2002l) for the bar fields would 
also add about 0.1 mag to the distance modulus. Figure |7| 
illustrates the effect of systematics in the assumed stellar 
populations: if we assume that the inner bulge is formed of 
an old population with enhanced a abundances, while the 
bar is formed of a younger population with solar abundances, 
the deduced distances show a much smaller dispersion from 
a galactic bar simple linear regression model. If confirmed by 
spectral observations, this implies a bar population which is 
formed from the inner disc, and not from the old bulge. 



4 THE STRUCTURE OF THE INNER 
GALAXY 

Figure [Ij] summarises our derived spatial map of the inner 
Galaxy. As well as the mean distances down each line of 
sight, the range of distances corresponding to the decon- 
volved one-cr dispersion in distance modulus is also illus- 
trated. In spite of this representation, we emphasise that 
the true distribution in distance has no reason to be Gaus- 
sian, as it is linked to the geometry of the bar/bulge and 
to our viewing angle. Our observation angle biases both our 
distance and dispersion estimates, because of the larger vol- 
ume, and hence probable larger number of stars present in 
the far side of our observation cone than in the near side. 
One should then keep in mind that our distance estimates 
will be slightly biased towards larger distances, especially 
for structures with a large true dispersion in distances, and 
as close to us as is observed at ^=+9°. 

Our mean red clump distances measured at £ = +9° and 
£ = —5° are consistent with the more recent bar models, par- 
ticularly those derived considering both gas dynamics and 
available surface brightness luminosity distributions. Those 
models deduce a bar orientation with respect to t he Sun - 
galac t ic centre directi on of 15° < (fib aT < 35° (e.g. [Ggrhard] 
l200ll lMerrifieldl 120031) . The model of iBissantz fc Gerhard! 
(2002) (their figure 11) is overlaid on figure 03 This model 
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Figure 6. The Galactic plane seen from the north Galactic pole. 
Solid circles indicate the mean positions of the red clump giants 
determined in this analysis at Galactic longitudes of i = ±9.7°, 
I = —5.7°, and t = 0°. Thick grey lines along the lines of 
sight through each mean distance represent the one— a range 
in distances deduced from a Gaussian fit to the red clump ap- 
parent distance modulus distribution, corrected for intrinsic red 
clump luminosity dispersion and photometric errors. The black 
line through the mean distances illustrates a Galactic bar of 3 kpc 
radius inclined at 22.5° to the Sun-Galactic Centre line. The back- 
ground contour map is a plane projection of the Bissantz & Ger- 
hard (2002) galactic bulge bar model, from their figure 11. The 
empty circles at I = +5.7° are only indicative, as the more distant 
point may be confused by background disk stars, as discussed in 
the text. The small cross at 3.5kpc from the Sun along the direc- 
tion i = +9.7° indicates the position of a local high-extinction 
region detected in our photometry. The square to the left of the 
dotted lines through our data indicates the mean position of the 
red clump detected by Hammersley et al. (2000) at ^=+27°. 

assumes a bar with an angle 4>b ar = 20°, a length of 3.5 kpc 
and an axis ratio 10:3-4. The overall agreement with our 
direct distance determinations is rather good. 

Our mean red clump distance derived at I — —9° is 
however not simply consistent with the bar model described 
above. A single linear fit to all the four mean distances mea- 
sured at I = +9, 0, —5, +9° has all four points more than 3<r 
from the 'best fit'. However the feature observed at £ = —9° 
does seem real: inspection of figure 0] shows a clear over- 
density of the whole red giant branch at the same distance 
as the mean clump distance, determined to be ~ 9 kpc, in 
this field. Moreover, in this line of sight, the most likely 
source of distance bias, the presence of disc red clump stars, 
would bias the distance estimate towards larger distances, 
not shorter. The red clump distance estimate at I = —9° is 
consistent with that at I = —5°. A possible physical inter- 
pretation is that we have detected the signature of the end 
of the bar. Star count peaks have also been detected in this 
region in studies of old OH /I R stars bv | Sevenster| ( 1999) and 
in the DENIS star counts bv | L6pez-Corredoira et al.l (2001). 
lLonez-Corredoira et al ] l|200ll) suggest that this local den- 
sity maximum may be related to the red clump overden- 
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Figure 7. Same as figure^ with the empty circles indicating the 
mean positions of the red clump giants determined with different 
isochrones for the bulge field at £ = 0° and the other fields. The 
solid circles distances assume that the stellar population of all the 
fields are 10 Gyr old and of solar metallicity. The empty circles 
positions are computed assuming for the bulge field at I = 0° a 10 
Gyr old population with enhanced a-element abundances, while 
assuming for all the other fields a 6 Gyr old population of solar 
abundances. 

sity detected by |HjMnme rslev et al] i2000T) at I = 27°. The 
iHammerslev et alT(|27)OoT ) density maximum is indicated in 
figure [S] by a square. 

The presence of these two local features at £ = 27° 
and £ = —9° and their mean red clump distance estimates 
could imply the presence of a double triaxiality in the inner 
Galaxy, a double bar, or a triaxial bulge oriented at ~ 22° 
with a longer, thiner bar oriented at ~ 44°. However, if such 
a second structure did exist, we should have detected its sig- 
nature in our other survey fields. No such complex signature 
is evident (cf. figure QJ. If two complex spatial distributions 
were projected down our lines of sight without being re- 
solved, this would imply that our distance determinations 
are biased so that the angle of the first structure is in fact 
smaller than the 22° we measured, and that the spread of the 
red clump distances observed down each line of sight would 
encompass both structures. The latter does not seem to cor- 
respond to the rather small distance dispersio ns we derived, 
which are summarised in table|5| Furthermore. IPicaud et al] 
(2003) do not detect a density excess at 1=-|-21 , confirm- 
ing that the structure seen at l=+27° is probably local, and 
unrelated to the larger bulge or bar. 

An alternative, and perhaps more consistent, interpre- 
tation of the £ = —9° structure we observe is the presence 
of a stellar ring or pseudo-ring at the end of the Galactic 
bar. If so, the bar and the stellar ring would have a radius 
of 2.3T0.25 kpc. Our observations at £ = +9° indicate that 
the radius of the bar is at least 2.7±0.2 kpc long. Those two 
determinations are consistent within one— a. We note that 
assuming the bulge/bar stellar population differences used 
for figure |7| leads to a larger bar radius of about 3 kpc. A 
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bar radius of about 2. 5 kp c would ag r ee wi th the model of 
ILepine fc Lerovl teOOOt) and ISevenstei lll999l) . but is smaller 
than the value indicated bv ToeTharc (j20fnf l in his review. 
The presence of a ring has been suspected by several au- 
thors. A molecular ring at about 4-5 kpc is a well known 
feature deduced from CO maps. Our detection of a step in- 
crease in the extinction distribution in field 9P (indicated 
by a cros s in fig ure El could be as sociated with this molec- 
ular ring. IComeron fc Torral |l996) link their derived distri- 
bution of ultracompact HII regions to the molecular ring, 
but also detect a star forming ring at about 2 kpc from the 
galactic centre. A stellar rin g was suspecte d to lie at 3.5 kpc 
from the galactic centre b ylBerte lli et alJ dl995T) . From an 
OH/IR star study ISevensterl " (^99ir "suggests that an inner 
ring lies between bet ween 2.2 and 3.5 kpc. From sta r counts 
in the DENIS survev lLopez-Corredoira et alJ i200ll) also ar- 
gue for the presence of a stellar ring, mainly from the detec- 
tion o f a density peak at I = —22° that, following ISevensterl 
( 1999), they associate with the tangential point to the 3- kpc 
arm, which is likely to be a (pseudo-)ring. The presence of 
a Galactic ring would also help to reproduce the obs erved 
microlensing optical depth iSevenster fc KalnaidEoOlT) . 

Inner ring structures are indeed frequently observed in 
barred spirals. Most inner rings are in fact pseudo-rings, 
formed from t he complex merging of inner bars and spiral 
structure (e.g. But a fc Combes! 1 1996f) . Such a pseudo-ring 
structure is consistent with all the various detections of sub- 
stru ctures between 2 and 4 kpc fro m the Galactic centre. 

lLopez-Corredoira et alJ (l200lT) detect an asymmetry in 
the Galactic plane extinction distribution, with more extinc- 
tion at negative than at positive longitudes. We do not con- 
firm this result. The fact that we do not find any asymmetry 
in red clump giant extinction between positive and negative 
longitudes (figure , while there is an asymmetry in dis- 
tance, is certainly consistent with a minimum in the dust 
distribution in this region. Interior of the bar outer radius a 
lower global density of stars and gas is often seen in observa- 
tions of other galaxies. Several authors have suggested that 
such a stellar density decrease may be present i n the Galactic 
disc teinnev et alJll99ll iBertelli et alJll995l. iFreudenreichl 
119981 ILepine fc Lerovl l200fj. lLopez-Corredoira et alJl2004l) 
Our observations are consistent with such a decrease in the 
dust distribution of the inner disc. 

Another interesting feature often observed in barred 
galaxies is the presence of dust lanes ass ociated with the bar . 
From DIRBE surface brightness maps, ICalbet et alJ <ll996t) 
suggested the presence of such a dust lane, preceding the 
Galactic bar. They estimate the extra absorption present 
at negative longitudes to be between 1 and 2.6 mag in K. 
The presence of such a feature is clearly ruled out by our 
observations. 

4.1 The line of sight towards £ = +5° 

The reddening-independent colour-magnitude diagram for 
the lines of sight towards I — +5° (figure |HJ differ signifi- 
cantly from those in other directions, in that the red clump 
stars indicate a very considerable range of distances, with no 
clear local maximum density. This is also very apparent in 
the relevant panels of figure which show the correspond- 
ing distribution in distance modulus. The broad apparent 
distance distribution observed for the red clump stars at 
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Figure 8. Colour - reddening independent magnitude diagram 
for field 5PP {£ = +5°). This figure is derived from H&K S pho- 
tometry, which have the faintest completeness limits for this field. 

I = +5° could be explained by the presence of two fea- 
tures, one at ~ 6 kpc and a second at ~ 11 kpc. Indeed, 
careful inspection of figure [8] suggests that the clump at 
K sH _ Ks ~ 12.5 does not extend redder than H-K s >1.2, 
while the fainter part, corresponding to distances greater 
than ~ lOkpc, extends beyond the photometric complete- 
ness limit. The first distance range is consistent with the 
bar structure confirmed abo ve. It is also consisten t with the 
fact that the OGLE data of IStanek et alJ Jl994), observed 
at a lower latitude of b=— 3.5°, show similar photometric 
behaviour at I = +5° compared to £ = —5°, except for the 
expected shift in distance modulus. 

The fainter feature in our low-latitude data could then 
be due to another structure, further away. This feature 
could be due to disc red clump stars. Indeed with in- 
creasing distance both the volume observed and the ex- 
tinction increase, and considering the logarithmic relation 
between distance and magnitude, any distant red clump 
stars visible would concentrate in the high-extinction high 
distance-modulus bottom-right part of the CMD. Deeper 
data than we have available are then needed to determine 
if we are indeed seeing the distant disk beyond the bulge 
or a local structure which could be associated with the 
ring. This region is also interesting for further study as a 
broad velocity-width m olecular clump has been detected at 
5.2° < £ < 6.0° (e.g. feovce fc Cohenl Il994l iBitran et alJ 
Il997t iDame : et alj|200ll), which m ay be a manifestation of 
gas shocks (iKumar fc Riffer3ll99llFuxlll999li . An estimate 
of the distance of those two structures, if the fainter one is 
real and indeed a local maximum, is presented in figureHJas 
the open circles. 



5 CONCLUSION 

We have obtained, reduced and analysed deep wide-field 
near-infrared photometry of five lines of sight towards the 
Galactic bulge within 0.25° of the Galactic plane, using 
the CIRSI camera. This has provided a substantial im- 
provement in our quantitative knowledge of inner Galactic 
structure, and in particular of the structure of the Galactic 
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bar. Our use of near-infrared data allows us to determine 
reddening-independent magnitudes and to observe bulge red 
clump stars within the galactic plane, where the extinction 
is too high for optical studies. This allows the detection of 
tracer stellar populations independently of their scale height. 
Colour - reddening-independent magnitude diagrams have 
been shown to disentangle the effects of distance and extinc- 
tion, allowing a direct conversion of red-clump star photom- 
etry (J — K,Kj-k) into an (Ay , distance) estimate. Those 
determinations have led to the following main results: 

- The presence of a triaxial structure at the centre of our 
Galaxy is confirmed. Its angle relative to the Sun-Galactic 
centre line is 4>bar = 22 ± 5.5° . It extends to at least 2.5 kpc 
from the Galactic Centre. A large axis ratio is excluded, but 
our data are consistent with a 10:3-4 ratio. In particular the 
distance dispersion of the bulge along the (£ =0°,b=l°) line 
of sight is less than 1 kpc. 

- A structure present at £ = —9.8° is not aligned with 
this triaxiality. We suggest that the structure present at £ = 
—9.8° is likely to be the signature of the end of the Galactic 
bar, which is therefore circumscribed by an inner pseudo- 
ring. 

- A decrease in the dust distribution inside the bar radius 
is inferred from the extinction distribution in our fields. 

- Our observations are not consistent with the existence 
of the dust lane prec eding the Galactic b ar at negative lon- 
gitudes suggested bv lCalbet et alJ il996f) . 
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